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SUMMARY  OF  RESULTS 


Significant  results  of  the  contract  study  are  summarized  below.  Detailed 
analytical  results  are  provided  in  Section  5,  and  illustrative  numerical 
results  are  given  in  Section  6. 

•  A  general  integral  expression  was  derived  relating  the  average 
efficiency  of  heterodyne  detection  to  the  degree  of  coherence  function, 
characterizing  the  spatial  phase  variations,  at  the  entrant  e  pupil  of  a 
generalized  laser  transceiver  system. 

•  Specific  integral  expressions  were  derived  for  the  heterodyne 
efficiency  of  a  transceiver  having  arbitrary  entrance-pupil  radius, 
receiver  focal  length,  local  oscillator  aperture,  and  local  oscillator 
Gaussian  amplitude  distribution  and  focal  length  in  response  to  a 
Schell-model  source  propagating  through  atmospheric  turbulence  and  an 
arbitrary  target  laser-reflectance  distribution  viewed  through 
atmospheric  turbulence. 

•  Computer  codes  were  developed  that  evaluate  the  laser  heterodyne 
efficiency  for  both  Schell-model  sources  and  general  target  reflectance 
functions. 

•  Computer  codes  were  developed  that  calculate  projected  laser  radar 
cross-section  distributions  for  both  spherical  and  conical  targets 
having  arbitrary  (homogeneous)  surface-roughness  characteristics. 

•  Heterodyne-detection  efficiencies  were  calculated  for  representative 
transceiver  and  target  parameters.  These  calculations  provide  specific 
functional  relations  for  the  Increase  in  heterodyne  efficiency  with 
increasing  target  range,  decreasing  target  size,  and  decreasing 
receiver  entrance  pupil  along  with  other  desired  functional  relations 
involving  other  transceiver  parameters  and  target  characteristics. 
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3.  CONCLUSIONS 

Salient  conclusions  resulting  from  this  study  are  summarized  below. 

•  Heterodyne-detection  inefficiency  resulting  from  target-induced  laser 
speckle  can  significantly  degrade  LADAR  performance  for  close  ranges 
(~  1  km)  and  fully  illuminated  target  conditions  unless  small  aperture 
receivers  (on  the  order  of  a  speckle  cell  diameter)  are  employed. 
Inefficiencies  of  IS  dB  or  greater  can  result  at  close  range. 

•  Atmospheric  scintillation  does  not  significantly  degrade  heterodyne- 
detection  efficiency  for  target  ranges  up  to  10  km  at  a  10.6  um  laser 
wavelength.  A  5  dB  degradation  in  slgnal-to-noise  ratio  is  the  maximum 
effect  expected  for  an  extended  target  at  10  km  range. 

•  Smooth  targets  generally  provide  a  higher  signal-to-noise  ratio  for 
heterodyne  detection  than  corresponding  rough  targets  because  of  the 
greater  spatial  coherence  at  the  receiver. 

•  A  uniform  local-oscillator  beam  (such  as  obtained  from  an  apertured 
TEMqq  laser  beam)  that  encloses  a  slightly  smaller  area  on  the  photo¬ 
detector  than  the  target  image  is  required  for  optimum  heterodyne 
detection  efficiency. 
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4.  RECOMMENDATIONS 

•  The  receiving  apertures  used  for  heterodyne  detection  of  laser  radar 
signals  should  be  designed  such  that  the  spatial  extent  of  the  aperture 
is  minimized  consonant  with  requirements  for  total  received  signal.  A 
circular  receiving  aperture  (non-Cassegrainian)  is  optimum  for 
minimizing  speckle-induced  heterodyne  inefficiencies  resulting  from 
laser  scattering  from  a  symmetric  diffuse  target. 

•  Laser  radar  transceivers  should  be  designed  so  that  the  local 
oscillator  beam  is  approximately  uniform  and  encloses  an  area  on  the 
photodetector  within  the  target  image. 

•  Experiments  should  be  performed  to  measure  the  efficiency  of  laser 
heterodyne  detection  for  a  variety  of  target  and  receiver 
configuration^  that  conform  to  the  analytical  assumptions  used  in  this 
study.  A  de  :ailed  comparison  between  measured  and  calculated  results 
should  be  performed  to  determir.  the  degree  of  validity  of  the  analyses 
and  any  inefficiency  effects  requiring  further  Investigation. 

•  The  anlytlcal  techniques  developed  in  this  study  should  be  extended  to 
treat  the  sum  and  difference  signals  arising  from  two  heterodyne 
detectors  in  the  focal  plane  of  a  receiver  system,  simulating  the 
tracking  signal  arising  from  a  raonopulse  LADAR.  Numerical  computations 
of  the  resultant  expressions  should  determine  the  magnitude  of  the 
tracking-error  signal  arising  from  speckle  and  scintillation  effects. 
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5.  ANALYSES 


Analytical  results  obtained  under  the  contract  study  are  presented 
below.  Section  5.1  provides  the  basic  heterodyne-efficiency  derivation,  while 
Sections  5.2  and  5.3  respectively  treat  the  efficiency  resulting  from  a 
Schell-model  source  and  a  generalized  laser  radar  cross-section 
distribution.  Section  5. A  summarizes  the  treatment  of  target  reflectivities. 

5.1  Laser  Heterodyne-Detection  Efficiency 

This  section  provides  an  analysis  of  the  average  efficiency  of  heterodyne 
detection  for  a  generalized  transceiver  system  in  response  to  a  phase  front 
that  is  distorted  by  diffuse  scattering  from  the  target  (speckle)  and 
atmospheric  turbulence.  Because  each  pulse  detected  by  the  transceiver 
represents  one  sample  of  statistically  equivalent  configurations  of  the  target 
and  the  atmosphere,  the  average  efficiency  is  treated  by  considering  the 
phase-front  correlations  (l.e.,  the  mutual  Intensity  function)  that  result 
from  an  average  over  the  ensemble  of  possible  configurations.  In  practice, 
wind  and  target  motions  provide  a  temporal  succession  of  pulse  samples  such 
that  a  time  average  closely  approximates  the  assumed  ensemble  average. 
Utilizing  the  development  of  Reference  1,  the  ensemble-averaged  photodetector 
current  correlation  function  of  a  heterodyne  LADAR  system  can  be  expressed  as 


<*IF(t)  xi F(t“T)) 


t*  IHh2' 


r(r,S,T) 
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heterodyne  intermediate-frequency  photocurrent 

scalar  mutual  coherence  function  incident  on  photodetector 

time 

transverse  coordinate  across  photodetector  surface 
local  oscillator  field 
pulse  shape 

Doppler-shifted  heterodyne  frequency 


5 


T  *  t  -  t' 

a  -  nGe/2hv 

n  «  detector  quantum  efficiency 

G  ■  gain 

e  *  electronic  charge 

h  *  Planck's  constant 

zQ  *  free-space  impedence 

v  *  photon  frequency. 

Equation  (1)  is  identical  to  Equation  (15)  of  Reference  1  except  that  one 
component  of  the  field  is  considered  (hence  the  scalar  form  of  the  general 
vector  expression  is  employed)  and  the  temporal  pulse  form  has  been  separated 
from  the  target/atmosphere-dependent  mutual  coherence  function.  Noting  that 
the  current  power  spectral  density  S(v)  and  the  current  correlation  function 
are  related  by  the  Fourier  transform  pairs. 


S(v) 


00 

J  dr  e~i2lTVT  <xip(t)  xip(t  -T))>, 


and 


00 

<XIF(t)  XIF(t  -  T)>  -  j 


,  i2nvx„  .  . 
dv  e  S  (v) , 


(2a) 


(2b) 


it  is  clear  that  the  ensemble-averaged  current  power,  S^,  integrated  over  the 
(relatively  small)  frequency  bandwidth  resulting  from  current  fluctuations  is 
given  by 


00 


~£0 


Thus,  the  frequency-integrated,  ensemble-averaged  current  power  depends  only 
on  the  current  variance.  Therefore,  the  temporal  variable  t  can  be  eliminated 
from  further  consideration  in  the  evaluation  of  with  the  understanding  that 
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The  mutual  coherence 


the  function  T  is  defined  only  when  the  pulse  is  "on." 
function  evaluated  for  zero  time  delay  (t  ■  0)  is  known  as  the  mutual 
intensity  function  (MIF)  and  will  be  referred  to  as  such  in  the  ensuing 
analysis. 

It  is  convenient  for  the  purposes  of  subsequent  analyses  to  define  the 
interference  integral  of  Equation  (1)  iu  a  plane  (having  transverse 
coordinate  f)  located  in  front  of  the  receiving  aperture.  Denoting 
by  (t)  the  local  oscillator  field  amplitude  back-propagated  from  the 
photodetector  to  the  receiving  aperture  plane,  the  average  frequency 
integrated  current  power  is 


nt,t  )  ”ib 


(t) 


(A) 


where  is  the  back-propagated  aperture  area  (i.e.,  the  minimum  areal 
subtense  of  either  the  receiving  aperture  or  the  back-propagated  local- 
oscillator  wave,  whichever  is  less).  Equivalently,  Equation  (A)  can  be 
written 


d^d2?  Y (|,|")  u£b(f)  u*b(|'). 


(5) 


where  u^b  is  the  normalized  spatial  variation  of  the  back-propagated  local 
oscillator  wave  and  U  is  the  magnitude,  viz., 


'  u*b  “«,<*>• 


(6) 


where  iu^COl^^  “  and  is  the  corresponding  normalized  MIF,  viz.. 


r (?,!')  -  uo2  y(I ,t'), 


(7) 
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where  |y(t,f")l  =  1.  When  the  average  intensity  [i.e.,  <U(f)  U*(f)>  •»  F(f,f')] 

max 

is  approximately  constant  over  the  receiver  aperture  (true  for  most  cases  of 
interest),  is  the  spatial  degree  of  coherence  function. 

Assuming  shot-noise-limited  detection,  the  frequency-integrated  current 
noise  power  is* 
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where  xD(,  is  the  direct  current  level  and  Av  is  the  bandwidth.  The  direct 
current  level  results  from  the  local  oscillator  power  level*,  and  hence  the 
noise  power  can  be  written  as 


NT  -  2eGAv  x  2a  P  «  -2-e— ■- 
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Thus,  the  average  current  power  signal-to-noise  ratio  is  given  by 
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Clearly,  optimal  signal-to-noise  ratios  are  achieved  when  the  incident  field 
replicates  the  spatial  distribution  of  the  local  oscillator  field,  i.e., 


-  u lh(t)  u t*d')  , 


yielding  the  ideal  signal-to-noise  ratio 


hvAv  * 


ideal 
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where 


|u£b(f')|2  (13) 

is  the  received  signal  power  for  a  perfectly  phase-matched  incident  signal  in 
the  receiving  aperture  encompassing  the  area  Ay*  In  general,  however,  the 
phase  fluctuations  occurring  in  the  incident  field  preclude  obtaining  the 
ideal  (maximum)  signal-to-noise  ratio  so  that  the  actual  ratio  achieved  as  a 
result  of  the  spatial  phase  fluctuations  can  be  expressed  as 
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.hi, 
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where  e  is  the  local  oscillator  efficiency  factor  defined  as 
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(15) 


Although  Equation  (15)  provides  a  valid  representation  of  the  local-oscillator 
efficiency,  it  is  often  more  appropriate  for  LADAR  analyses  to  consider  the 
maximum  signal-to-noise  ratio  that  is  achievable  from  the  total  power  enclosed 
by  the  receiving  aperture,  i.e., 


max 
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where 
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(16) 


(17) 
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and  it  is  again  assumed  that  the  received  intensity  is  constant  over  the 
receiving  aperture.  In  this  case  the  maximum  power  signal-to-noise  ratio 
heterodyne  efficiency  factor  is  given  by 


^ - — — - — j -  .  (18) 

\  ffdt  ''V^1 

V 

5 . 2  Schell-Model  Source  Efficiency  Factors 

To  permit  specific  calculations  of  LADAR  efficiency  factors,  the 
generalized  transceiver  system  depicted  in  Figure  1  is  considered  and  a 
Schell-model  source  is  assumed  for  the  ensuing  development.  The  parameters  of 
the  Figure  1  transceiver  are  as  follows? 

A  -  local  oscillator  aperture  radius, 

P  «  entrance  pupil  radius, 

-  focal  length  of  local  oscillator  lens/mirror  (lens  3),  and 
f r  *  focal  length  of  receiving  lens/mirror  (lens  2) . 


10 


i  TEMqq  mode  local  oscillator  beam  is  assumed  such  that  the  (scalar) 
local  oscillator  amplitude  is  given  by 


VS)  ‘  'V  e'li5|2/L“ 


in  the  plane  (having  transverse  coordinates  p)  located  in  front  of  lens  3. 

Because  a  Schell-model  source  is  assumed,  the  incident  mutual  intensity 
function  can  be  expressed  as 


.  ,  (a|||2-ta*|r|2)  -|l4'I2/r  2 

r  (!,!')  -U/e  1  e  c 


where  rc  is  the  coherence  radius  and 


„  .  2  ,  i 

*  ^  +  V 


where  LA  is  the  amplitude  radius  and  is  the  average  phase-front  radius  of 
curvature. 

Utilizing  the  Fresnel  diffraction  integral,  the  local  oscillator  field 
amplitude  in  the  photo-detector  plane  is  given  by 


V*r) 
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where 


+  2  2 
1  when  p  <  A 

*►2  2 
0  when  p  >  A 
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is  the  local-oscillator-pupil  function.  Similarly,  the  photo-detector  field 
is  related  to  the  back-propagated  local-oscillator  field  by  the  expression 
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Hence,  the  back-propagated  field  can  be  expressed  as 
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where 


ikr2 


H(r)  -  e 


(26) 


and 


(27) 


utilizing  the  inverse  Fourier  transform  of  Equation  (24).  Equation  (25)  can 
be  expressed  in  convolution  form,  i.e., 
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d2t'  H  (t  -  I')  G(f')  , 


(28) 


12 


where 


(29) 


Thus,  the  back-propagated  local-oscillator  field  is  given  by 


(32) 

where  0(x)  is  the  Heaviside  function. 

The  Integral  form  of  Equation  (32)  is  not  analytically  integrable  because 
of  the  finite  range  imposed  by  the  Heaviside  function.  If  it  is  assumed  that 
the  aperture  area  goes  to  infinity.  Equation  (32)  is  easily  evaluated, 
yielding 
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(33) 


(34) 


(35) 


A  result  that  is  identical  to  Equation  (33)  is  obtained  if  Equation  (32)  is 
evaluated  by  the  method  of  double  stationary  phase, ^  providing  the  stationary 
phase  point,  given  by 


l,  m  |  (1  -  4iz"L'**'2) 

’°  1  +  (4z'I/“2)2  * 


(36) 


lies  within  the  domair  of  integration.  Because  of  the  high  frequencies 

-  '-2 

associated  with  LADA l  systems,  the  factor  z  L  appearing  in  Equation  (36)  is 

much  less  than  unity  (it  is  zero  for  matched  focal  lengths,  i.e.,  f.  m  f  )  and 

x*  r 

thus  %'  "  *•  Because  the  stationary  phase  result  equals  the  infinite  limit 
result  (Equation  (33)],  the  major  contribution  to  the  Integral  arises  from  the 
domain  of  integration  in  the  immediate  neighborhood  of  the  stationary  phase 
point.  Therefore,  the  integral  of  Equation  (32)  is  approximately  given  by 


Utilizing  the  Equation  (37)  result  for  the  back-propagated  local-oscillator 
wave,  the  signal  power,  Equation  (5),  can  be  expressed  as 
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where 
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and  Che  Schell-model  source  of  Equation  (20)  Is  assumed.  Equation  (38)  can  be 
written  in  more  compact  form,  viz., 
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utilizing  the  properties  of  the  Heaviside  function.  Equation  (42)  can  now  be 
approximately  integrated  utilizing  the  double  stationary-phase  method  yielding 
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The  parameters  of  many  problems  of  interest,  however,  invalidate  the 
stationary-phase  result  of  Equation  (45)  (see  Section  6),  so  that  in  many 
cases  numerical  integration  techniques  must  be  employed.  A  transformation  to 
polar  coordinates  permits  analytical  integration  over  the  angular  variables  so 
that  Equation  (42)  can  be  reduced  to  a  two-fold  integral  of  the  form 


ST  -  U  2  | K | 2  (2n)2 
o 


R_  R 
m  m 

// 

*^o  o 


dcdrcr  i 


.2  2 
-  a  K  -a*  K 
o  o 


(46) 


where  IQ  is  the  zeroth-order  modified  Bessel  function  of  the  first  kind.  The 
finite  Integrals  of  Equation  (46)  can  be  integrated  by  standard  numerical 
methods.  The  Integral  result  of  Equation  (46)  can  be  combined  with  the 
relatively  straightforward  integration  of  the  required  local  oscillator 
Intensities,  using  the  back-propagated  wave  result  [Equation  (37)],  to  obtain 
the  following  specific  functional  forms  for  the  local-oscillator  efficiency 
[Equation  (15)]  and  maximum  heterodyne  efficiency  [Equation  (18)]: 
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5 . 3  Heterodyne  Detection  Efficiency  for  General  LRCS  Distributions 


The  vacuum,  far-field,  mutual  intensity  function  resulting  from 
Incoherent  laser  scattering  from  a  rough  object  can  be  written^ 


_  lit  t+  ±- 

c  (T.  n  ~ . *  • 

where 

t*  -  {  (I  +  !')  , 

r  -  ?  -  r. 

R  -  range  to  target  , 
p  •  transverse  coordinate  to  line-of-sight  at  the  target  plane  , 

and 

2  •*  ^ 

■  differential  projected  scattering  cross-section  per  unit  area. 
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Equation  (49)  replicates  Equation  (1)  of  Reference  1  except  for  constant 
factors,  the  conversion  of  angle  coordinates  to  linear  displacements,  and  the 
rewriting  of  the  range  difference  factor  in  terms  of  sum  and  difference 


coordinates  (i.e.,  R  -  R"  ■  i+*t  /R) •  Although  the  Equation  (49)  prescription 
for  the  M1F  is  valid  only  in  the  far-field,  the  Reference  5  development 
demonstrates  that  this  range  is  much  smaller  than  the  classic  far-field  range 
(on  the  order  of  the  body  dimensionality)  for  relatively  incoherent 
scatterers.  Reference  6  demonstrates  that  atmospheric  turbulence  modifies  the 
vacuum  MIF  in  the  following  manner: 
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where  5  and  s  are  unit  pointing  vectors  and  F  is  the  Fourier-transformed  two- 
source,  spherical-wave,  atmospheric  MIF.  Because  the  integral  of  Equation 

(49)  is  functionally  independent  of  the  average  field-point 

coordinate,  the  integrated  functional  form  for  the  vacutra  MIF  in  Equation 

(50)  is  that  of  the  exponential  multiplying  the  integral  of  Equation  (49), 
which  is  the  functional  description  of  a  spherical  phase  front.  Thus,  the  MIF 
that  describes  the  field  correlations  resulting  from  incoherently  scattered 
laser  radiation  propagating  through  atmospheric  turbulence  is  of  the  form 
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2  2-3/5  ft 

where  p  »  (3/8  x  1.455  k  RC  ]  '  is  the  spherical  wave  coherence  length. 

0  n 

Glint  points  on  the  surface  of  an  object  scatter  radiation  in  the  form  of 
spherical  waves  with  a  scattering  cross  section  (in  the  geometrical  optics 
approximation)  given  by  it  times  the  product  of  the  Gaussian  radii  of 
curvature.  Thus,  assuming  (for  the  purposes  of  subsequent  analyses)  a  rough 
object  having  one  centrally  located  glint  point,  the  degree  of  coherence 
function  will  be  given  by 
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where  a^  and  a£  are  the  Gaussian  radii  of  curvature  at  the  glint  point,  A  is  a 

factor  that  relates  the  relative  strength  of  the  coherent  (glint)  and 

incoherent  components,  and  the  denominator  of  Equation  (52)  results  from  the 

normalization  requirement  of  the  degree-of-coherence  function.  Reference  7 

demonstrates  that  the  coherent  component  of  surface  scattering  is  attenuated 

by  an  exponential  factor  of  the  form  where  represents  the  phase 

variance  and  can  be  related  to  the  rough-surface  height  deviation.  Reference 

8  shows  that  rough-surface  speckle  statistics  are  replicated  when  the 

incoherent  component  of  the  MIF  is  attenuated  as  the  phase  deviation  decreases 

2 

according  to  the  factor  (1  -  e”0^  ).  Thus,  the  relative  strength  factor  A  can 
be  estimated  from  the  phase  variance  at  the  target  according  to  the  formula 


A  ■ 


(53) 


which  is  also  consistent  with  the  MIF  of  Reference  2. 

Recognizing  that  the  generalized  degree  of  coherence  function  of  Equation 
(52)  consists  of  the  sum  of  a  coherent  and  incoherent  component,  the 
calculated  heterodyne  efficiency  factor  given  by  Equation  (18)  can  likewise  be 
split  into  the  sum  of  two  components.  The  maximum  heterodyne  efficiency  can 
be  expressed  as 
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where  o  is  the  total  target  incoherent  LRCS,  i.e., 
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and  the  coherent  and  incoherent  efficiencies  are  determined  from  the  relation 
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with  the  integral  factors,  Qcoh/incoh>  <*e*ined  by 
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and 


where  the  subscripts  r  and  i  on  the  parameter  8  denote  the  real  and  imaginary 
parts,  respectively.  Because  the  coherent  component  of  the  degree-of- 
coherence  function  is  in  the  form  of  a  Schell-model  MIF,  the  results  of 
Section  5.2  can  be  used  to  express  the  coherent  Integral  factor  in  the  more 
condensed  form 


To  numerically  evaluate  the  incoherent  component  of  the  heterodyne 
efficiency,  it  is  necessary  to  simplify  the  incoherent  integral  factor, 
Equation  (58).  Noting  the  Fourier  kernel  in  the  integral  of  Equation  (58), 
the  convolution  theorem  can  be  employed  to  write 
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where  the  double  Fourier  transform  properties  of  a  function  have  been  utilised 
to  write  Equation  (64).  Employing  a  change  of  variables 
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so  that 
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Equation  (61)  is  written  in  the  form 
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where  the  function  F  filters  the  spatial  frequencies  of  the  power  spectral 

density  of  field  fluctuations^  represented  by  the  scattering 
2  +  +2  »v 

distribution  d  o(x)/dx  •  The  filter  function  F  describes  the  limited  bandpass 
of  spatial  frequencies  passed  by  the  atmosphere,  receiving  and  local 
oscillator  apertures,  and  the  local  oscillator  amplitude  distribution. 

Although  the  convolution  theorem  can  be  applied  to  the  integral  of  Equation 
(63)  to  write  the  filter  function  in  one-dimensional  form  (after  an 
Integration  over  an  angle  variable)  as 
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it  is  desireable  to  obtain  an  approximate  analytical  formula.  Appendices  A 
and  3  provide  derivations  for  the  following  approximate  formulae: 
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Numerical  calculation  (illustrated  in  Section  6)  of  the  filter  function  and 
the  approximations  given  in  Equations  (70)  and  (71)  show  that  the  following 
combination  of  these  approximations  provides  an  excellent  estimate  of  the 
exact  expression  [Equation  (68)]: 
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It  is  easy  to  demonstrate  using  the  approximate  forms  for  the  filter  function 
that 
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as  required  from  the  physics  of  the  problem. 
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5.4  Target  Scattering  Properties 

5.4.1  Spherical  Target  Schell-Model  Approximation 

Because  of  the  approximate  analytical  results  [Equation  (45)3  obtained 
for  the  heterodyne  signal  resulting  from  a  Schell-model  MIF  and  the  relative 
simplicity  of  the  exact  Schell-model  efficiency  expressions  [Equations  (47)  - 
(48)],  it  is  desireable  to  obtain  an  expression  for  the  Schell-model  MIF  at 
the  receiver  aperture  that  simulates  a  target  scattering  problem  of 
interest.  Spheres  provide  excellent  targets  for  laser  radar  experiments 
because  the  scattering  is  independent  of  target  orientation.  The  following 
development  provides  expressions  for  a  Schell-model  MIF  that  approximates  the 
MIF  arising  from  scattering  from  spherical  targets. 

Reference  2  provides  a  derivation  of  the  MIF  arising  from  a  partially 
coherent  beam-wave  source  propagating  through  atmospheric  turbulence.  The 
derived  field-point  MIF  is  given  by 
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where  and  are  the  y  components  of  the  sum  and  difference  field-point 
coordinates,  i.e., 
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and 
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so  that 
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where 
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The  MIF  is  calculated  in  terms  of  the  y  component  (i.e.,  •  $2x  -  0) 

because  it  contains  the  interesting  angular  (9)  dependence.  The  effective 

A 

amplitude  (L')  and  curvature  radii  (L^)  as  well  as  the  coherence  lengths 
(p')  in  Equation  (75)  are  given  by: 
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F  -  1  -  e  *  ,  (88) 
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D2  -  { [p~2  Z2  +  |1  +  iaR | 2 ]  [z2p“2  cos20  +  |1  +  iaR|2]}  .  (90) 

Approximate  values  of  the  source  amplitude  radius  (L^)  and  phase  curvature 

radius  (L  )  are  derived  below  such  that  the  asoumed  Schell-model  source 
9 

emulates  the  known  scattering  properties  of  a  sphere. 

The  scattering  cross  section  of  a  target  is  defined  by  the  relation 
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The  scattering  cross  section  of  a  perfectly  specular  sphere  is  given  by 
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(where  Lgp  is  the  sphere  radius)  when  the  sphere  is  large  with  respect  to  a 
wavelength.  Thus,  a  necessary  requirement  for  the  Schell-model  source  sphere 
simulation  is 
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The  sphere  radius  has  been  used  for  the  phase  front  curvature  in  Equation  (93)  1 
because  of  the  known  curvature  of  a  wave  front  scattered  by  a  sphere.  The  j 
Equation  (93)  requirement  yields  j 


(94) 


The  predicted  amplitude  radius  of  Equation  (94)  is  approximately  equal  to  the 
first  Fresnel  zone  of  the  sphere  which  is  consistent  with  the  physical  optics 
prediction  of  specular  scattering  from  a  sphere. 

The  cross  section  associated  with  the  incoherent  (i.e.,  diffuse) 
component  of  scattering  from  a  sphere  can  be  estimated  by  integrating  the 
diffuse  differential  cross  section  over  the  observed  sphere  surface  area. 
Because  the  cross  section  per  unit  illuminated  area  (o°)  is  related  to  the 
bidirectional  reflectance  distribution  function  (p'  )  by  the  relation'* 
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the  diffuse  sphere  cross  section  is 
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If  a  surface  scatters  all  the  energy  incident  upon  it  in  accordance  with 
Lambert's  law  with  half  of  the  scattered  power  in  each  of  two  orthogonal 
polarization  components,  the  bidirectional  reflectance  distribution  function 
is  a  constant  equal  to  (2rr)"”*.  Thus,  the  cross  section  associated  with  a 
Lambertian  spherical  scattering  surface  is 


bpRefetian 


sin0  cos2  8d8  -  4  2 

3  sp 


Because  the  Lambertian  sphere  cross  section  is  within  a  factor  of  1/3  of  the 

specular  sphere  cross  section,  it  Is  reasonable  to  estimate  the  cross  section 

2 

for  the  diffuse  scattering  cross  section  by  regardless  of  the  reflectance 

distribution  function  of  the  surface.  Thus  a  Schell-model  approximation  of 
the  scattering  MIF  resulting  from  a  spherical  target  is  obtained  by  using  the 
relations 
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in  Equation  (75)  along  with  the  requirements  that 
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in  Equations  (81)  -  (87).  The  final  relationship  [Equation  (99c)]  arises  from 
the  requirement  that  the  amplitude  radius  for  the  diffuse  scattering  process 
must  be  the  the  same  order  of  magnitude  as  the  actual  sphere  radius. 

A  check  on  the  validity  of  the  spherical  source  Schell-model 

approximation  is  obtained  by  calculating  the  coherence  length  arising  from  the 

2 

incoherent  component  without  atmospheric  turbulence  (i.e.,  »  0) .  It  is 

straightforward  to  show  that  for  a  range  much  larger  than  a  sphere  radius  and 
a  target  coherence  length  less  than  a  wavelength. 
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This  result  can  be  compared  with  the  coherence  length  arising  from  a  uniformly 
illuminated  disc  (corresponding  to  the  brightness  distribution  of  a  Lambertian 
sphere).  The  degree-of-coherence  function  for  a  uniform  disc  source  is 
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The  e“*  value  of  the  Equation  (102)  distribution  occurs  for  an  argument  of  the 
Bessel  function  of  rpproximately  2.5.  Thus  the  Lambertian  sphere  coherence 
length  of 
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is  in  good  agreement  with  the  Schell-model  approximation  result.  It  is  also 
straightforward  to  show  that  the  predicted  coherence  length  of  the  coherent 
component  correctly  models  coherent  spherical  wave  propagation  through  the 
atmosphere. 

5.4.2  Projected  LRCS  Distributions  for  Spheres  and  Cones 

The  calculation  of  the  heterodyne  efficiency  resulting  from  Incoherent 
target  scattering  [Equation  (67)]  requires  the  functional  description  of  the 
projected  LRCS  distribution.  If  the  laser  bidirectional  reflectance 
distribution  function  (or  equivalently  the  cross  section  per  unit  area)  of  the 
target  material  is  specified,  the  LRCS  distribution  can  be  determined  from  the 
target's  geometry.  The  development  below  illustrates  the  calculation  of  LRCS 
distributions  for  two  simple  target  geometries:  spheres  and  cones.  Subsequent 
numerical  calculations  (presented  in  Section  6)  utilize  these  expressions 
together  with  the  ROSSCO  laser  scattering  algorithm^  to  specify  LRCS 
distributions  for  illustrative  numerical  calculations. 


Spheres:  The  projected  differential  area,  dA^,  resulting  from  a  surface 
differential  area,  dA,  that  is  tilted  an  angle  0  from  normal  with  respect  to 
the  projection  direction  is 


dA  “  cos0dA. 
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(104) 


Because  a  differential  surface  element  of  a  sphere  that  is  a 

distance  |?|  radially  outward  from  the  projection  normal  (see  Figure  2)  is 

tilted  at  an  angle  sin”*  (|?|/L  ),  the  distributed  LRCS  is  defined  by  the 
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where  the  argument  of  the  cross  section  per  unit  illuminated  area,  o°,  in 
Equation  (105)  is  the  backscatter  angle. 
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Figure  2.  Sphere  projection  geometry. 


Cones:  To  determine  the  LRCS  distribution  of  a  conical  surface,  the 
projection  geometry  depicted  in  Figure  3  is  employed  where  Xp  and  z p  denote 
the  orthogonal  coordinates  of  the  projection  plane.  The  surface  of  the  cone 
is  defined  by  the  eqv.ation 


(106) 


while  the  polar  angle,  u,  of  a  point  on  the  surface  of  the  cone  is  given  by 
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Figure  3.  Cone  projection  geometry. 


Because  the  projection  coordinate  Zp  is  related  to  the  cone  surface 
coordinates  by  the  transformation  relation 


z  ■  z  sin0  +  y  cosB, 
P  1 


where  0  is  the  projection  direction,  the  equation 
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relates  the  projection  plane  coordinates  to  the  y  body  coordinate  on  the 
surface  of  the  cone*  Equation  (109)  has  two  solutions  for  a  given  projection 
plane  coordinate  (Xp,  Zp)  corresponding  to  the  two  intersection  points,  on  the 
near  and  far  side,  of  the  infinite  cone.  The  scalar  product  of  the  projection 
direction  vector,  n^,  with  the  cone  surface  nortral  direction  vector  yields  the 
following  equation  for  the  local  surface  tilt  angle,  0^ : 


cos0£  «  sin@  cosw  cosifi  +  cosB  sin<f>. 
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where 
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Thus,  the  projected  LRCS  distribution  is  given  by 
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where  the  (x*,  z*)  coordinates  fall  within  the  projected  rectangular  outline 
of  the  cone  defined  by  the  equations 
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the  projection  point  lies  on  the  surface  of  the  truncated  cone  defined  by 
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and  the  y(x^,  z )  solution  of  Equation  (109)  used  to  solve  Equations  (107) 
and  (110)  Is  the  minimum  of  the  two  roots. 
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6.  NUMERICAL  RESULTS 


To  establish  some  of  the  functional  relations  of  the  expressions  derived 
in  Section  5,  the  results  of  several  numerical  computations  are  presented 
below  in  graphical  and  tabular  form.  These  results  illustrate  the  effects  of 
various  transceiver  and  target  scattering  parameters  on  predicted  heterodyne 
performance.  The  various  parameters  chosen  for  these  illustrative 
calculations  were  selected  to  provide  representative  results  rather  than  to 
model  a  specific  transceiver  configuration  or  target  scattering  scenario.  A 
wavelength  of  10.6  ym  was  assumed  for  all  computations.  Because  calculated 
results  were  computed  with  18  different  computer  codes  that  were  developed  for 
specific  computations  in  the  course  of  an  evolving  analytical  research 
program,  the  codes  are  not  listed  in  this  report. 

6.1  ileterodyne  Signal  Calculations 

Figures  4-8  illustrate,  in  three-dimensional  graphical  form,  the 

variation  in  the  relative  heterodyne  signal  predicted  by  the  stationary-phase 

approximation  result.  Equation  (45),  (Figures  4a-8a) ,  and  the  exact  two- 

dimensional  integral  result.  Equation  (46)  (Figures  4b-8b).  Although  Figures 

4  and  5  demonstrate  surprisingly  close  agreement  between  the  stationary-phase 

result  and  the  exact  result  as  a  function  of  the  local  oscillator  parameters  A 

and  Lq  and  the  coherence  length,  rc,  only  general  features  are  replicated  in  the 

two  calculations  in  Figures  6  and  8  (for  the  variables  f . ,  L, ,  and  r  ),  and  a 

9  c 

totally  different  functional  relation  between  the  focal  lengths  f^  and  f^  is 
shown  in  Figure  7.  Thus,  although  the  stationary-phase  result  [Equation  (45)] 
provides  an  order-of-magnitude  estimate  for  the  signal  power  in  an  analytic 
form  which  can  be  quickly  calculated  on  a  hand  calculator,  it  cannot  be  used 
to  provide  accurate  predictions  or  uniformly  valid  functional  relations. 

The  general  functional  trends  depicted  in  the  exact-result  signal  power 
calculation  of  figures  4b-8b  can  be  interpreted  from  the  known  physics  of  the 
detection  process.  Figures  4-6  show  that  an  increased  coherence  length,  rc, 
of  the  detected  signal  always  results  in  an  increased  signal  power  since  the 
signal  becomes  more  coherent  (the  speckle  cell  dimensionality  increases)  and 
hence  the  phase  front  becomes  more  uniform  and  amenable  to  the  heterodyne 
detection  process.  Figures  4  and  5  also  demonstrate  the  increase  in  signal 
power  as  the  local-oscillator  amplitude  radius  or  aperture  radius  increases. 
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reflecting  the  increase  in  the  photodetec*"or  area  used  for  the  detection 
process.  The  oscillatory  behavior  shown  in  Figure  5  as  a  function  of  the 
local  oscillator  radius  (for  a  coherent  wave)  replicates  the  "ringing" 
interference  behavior  predicted  by  earlier  investigations.^’ ^  Figures  6b-8b 
shown  an  increase  in  the  predicted  heterodyne  signal  with  decreasing  local- 
oscillator  focal  length  because  the  receiving  aperture  spot-size  is  larger 
than  the  local-oscillator  spot-size  on  the  photodetector  (for  the  assumed 
local-oscillator  radius) ,  resulting  in  incomplete  heterodyne  detection  of  the 
signal  enclosed  by  the  receiving  aperture.  Equivalently,  the  analysis  of 
Section  5.2  shows  that  the  back-propagated  wave  underfills  the  receiving 
aperture  when  the  local-oscillator  radius  is  smaller  than  required  for 
complete  heterodyne  detection  of  the  received  signal.  It  is  noteworthy  that 
Figure  8  shows  that  the  heterodyne  signal  is  relatively  unaffected  by  the 
phase-front  curvature  of  the  received  signal. 


(b) 


Figure  4.  Later  heterodyne  signal  power  results  calculated  by  the  stationary  phase 

method  (a)  and  numerical  integration  technique  (b)  as  a  function  of  the  local 
oscillator  amplitude  radius,  Lq,  and  Schell-model  coherence  length,  rc.  Fixed 
parameters  of  the  calculation  are  A  *  10  cm,  P  *  9.5  cm,  fg  m  12.5  cm, 
fr  *  25  cm,  “  -200  m,  and  R  *  200  m. 
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Figure  5.  Later  heterodyne  signal  power  results  calculated  by  the  stationary  phata 
method  (a)  and  numerical  integration  technique  (b)  as  a  function  of  the 
local  oscillator  aperture  radius.  A,  and  Schell-model  coherence  length,  rc. 

Fixed  parameters  of  the  calculation  are  1.0  ■  8  cm,  P  *  10  cm,  fg  ■  12.5  cm, 
f  r  *  25  cm,  L^  ■  -200  m.  and  R  ■  200  m. 
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Figure  6.  Laser  heterodyne  signal  power  results  calculated  by  the  stationary  phase  (a) 
and  numerical  integration  technique  (b)  «s  a  function  of  the  local 
oscillator  focal  length,  fg,  and  Schell-model  coherence  length,  rc.  Fixed 
parameters  of  the  calculation  are  A  ■  2  mm,  P  *  2.5  cm,  L0  *  8  cm, 
fr  =  25  cm,  L ^  =  -200  m,  and  R  =  200  m. 


40 


OPOJ4WM-7 


Figure  8.  Laser  heterodyne  signal  power  rasults  caiculatad  by  tha  stationary  phase  (a) 
and  numsricaf  integration  technique  (b)  as  a  function  of  tha  local 
oscillator  focal  length,  fg,  and  phase  front  curvature,  L^.  Fixed  parameters 
of  tha  calculation  are  A  *  2  mm,  P  «  2.6  cm,  L0  *  2.5  cm,  fr  •  25  cm, 
rc  »  1  cm,  and  R  *  200  mm. 


6.2  Heterodyne  Efficiency  for  Schell-Model  Sources 


Figures  9-11  show  representative  results  for  the  calculated  maximum 
heterodyne  efficiency.  Equation  (48),  as  a  function  of  the  local-oscillator 
parameters,  A  and  LQ,  the  coherence  length,  rc,  and  the  focal 
length  f^  and  f^.  Similar  results  are  depicted  in  Figures  9  and  10: 
heterodyne  efficiency  decreases  with  decreasing  effective  local-oscillator 
aperture  and  coherence  length.  Figures  9  and  10  show  that  increasing  the 
local-oscillator  amplitude  radius,  LQ,  (for  sufficiently  large  aperture)  or 
increasing  the  local-oscillator  aperture,  A,  (for  sufficiently  large  amplitude 
radius)  will  result  in  an  increased  detection  efficiency  until  a  condition  is 
obtained  where  the  back-propagated  local-oscillator  beam  overfills  the 
receiving  aperture,  after  which  no  efficiency  benefits  are  obtained.  The 
primary  difference  between  Figures  9  and  10  is  the  more  gradual  increase  in 
the  heterodyne  efficiency  with  increasing  local-oscillator  amplitude  radius 
vis-a-vis  aperture  size.  The  Figure  11  plot  of  the  heterodyne  efficiency  as  a 
function  of  the  local-oscillator  and  receiver  focal  length  again  illustrates 
the  requirement  that  the  back-propagated  local-oscillator  beam  must  overfill 
the  receiving  aperture  for  optimum  signal-to-noise  ratio  efficiency.  Thus 
optimum  heterodyne  efficiency  is  generally  obtained  when  the  local-oscillator 
amplitude  radius  is  large  and  the  local-oscillator  aperture  truncates  the 
Gaussian  beam  profile  such  that  the  image  of  the  aperture  at  the  receiving 
lens  is  larger  than  the  entrance  pupil. 

Figure  12  provides  a  graph  of  the  computed  maximum  heterodyne  efficiency, 
Equation  (48),  as  a  function  of  the  sphere  radius  and  source  coherence  length 
used  in  the  Schell-model  approximation  of  the  radiation  characteristics 
arising  from  laser  scattering  by  spheres  (Section  5.4.1).  Whereas  previous 
Figures  (4,  5,  6,  9,  10)  show  the  decrease  in  heterodyne  efficiency  with 
decreasing  coherence  length,  Figure  12  depicts  the  coherence  length  effects  on 
the  efficiency  resulting  from  target  parameters.  The  source  coherence 
length,  p^,  used  in  Schell-model  sphere  approximation  parameterizes  the 
coherence  length  arising  from  the  laser  scattering  process  and  thus  relates  to 
the  surface-height  correlation  length  for  a  rough-surface  scattering 
process.  Clearly,  large  source  coherence  lengths  (much  larger  than  a 
wavelength)  describe  scattering  from  smooth  surfaces,  whereas  small  source 
coherence  lengths  result  from  diffuse  (rough  surface)  scattering  processes. 
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The  results  plotted  in  Figure  12  demonstrete  that  maximum  heterodyne 
efficiency  is  achieved  for  sphere  scattering  when  either  the  sphere  radius  is 
small  or  the  surface  is  very  smooth.  It  is  also  noteworthy  that  Figure  12 
indicates  that  large  spheres  must  be  smoother  than  smaller  spheres  to  achieve 
comparable  heterodyne  efficiency. 
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Figure  9.  Calculated  later  heterodyne  maximum  detection  efficiency  at  a  function  of 
the  local  oscillator  aperture  radius.  A,  and  Schell-model  coherence  length,  rc. 
Fixed  parameters  of  the  calculation  are  L0  »  1.25  cm,  P  *  1.25  cm, 
fj  ■  12.5  cm,  ff  “  25  cm,  “  -200  m,  and  R  -  200  m. 
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Figure  10.  Calculated  laser  heterodyne  maximum  detection  efficiency  as  a  function  of 
the  local  oscillator  amplitude  radius,  L0,  and  Schell-model  coherence 
length,  rc.  Fixed  parameters  of  the  calculation  are  A  »  1.25  cm, 

P  =  1.25  cm,  ffi  «  12.5  cm,  fr «  25  cm,  i.0  »  -200  m,  and  R  -  200  m. 
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Figure  11.  Calculated  laser  heterodyne  maximum  detection  efficiency  as  a  function  of 
the  local  oscillator  focal  length,  fg,  and  receiver  focal  length,  fr.  Fixed 
parameters  of  the  calculation  are  L0  »  1 .25  cm,  P  *  1.25  cm,  A  *  1.25  cm, 
rc  »  1  m,  =  -200  m,  and  R  «  200  m. 


QP03-OSSS-S 

Figure  12.  Calculated  laser  heterodyne  maximum  detection  efficiency  as  a  function  of 
a  spherical  targets  radius,  L$p,  and  scattered  laser  coherence  length,  p(J.  Fixed 
parameters  of  the  calculation  are  P  =  2.5  cm,  Lo  =  2.5  cm,  A  “  2.5  cm, 
ff  *  25  cm,  fg  =  25  cm,  H  =  500  m,  C„2  =  lo-1®  m'2/3  =  5 ,0.  The 

illuminating  later  beam  is  assumed  to  have  a  focal  length  of  -25.0  m  and  a 
2.5  cm  aperture  radius. 


j 
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6.3  Heterodyne  Efficiency  for  General  LRCS  Distributions 


6.3.1  Filter  F motion  Approximation 

To  account  for  the  effects  of  atmospheric  turbulence  and  finite  local 
oscillator  amplitude  radius  in  the  general  calculation  of  the  heterodyne 
efficiency  [Equations  (56)  and  (58)],  it  was  necessary  to  develop  an 
approximate  form  for  the  filter  function  [Equation  (72)]  to  reduce  the  exact 
six-fold  dimensional  integral,  Equation  (58),  to  a  more  numerically  tractable 
four-fold  form  [Equation  (67)].  Figures  13  and  14  provide  graphs  of  (a)  the 
exact  filter  function  [Equation  (68)3 ,  (b)  the  stationary-phase  approximate 
result  [Equation  (70)]  (Appendix  A),  and  (c)  the  series-expansion  approximate 
result  [Equation  (71)]  (Appendix  B)  as  a  function  of  the  unitless  variable 

r\ 

p  ■  R  )j'  and  the  scaling  parameter  8  **  R  for  a  relatively  short  range 

m  m 

of  p  (Figure  13)  and  a  more  extended  range  of  p  (Figure  14) .  The  results 
plotted  in  Figures  13  and  14  demonstrate  that  the  series-expansion  result 
replicates  the  exact  result  for  values  of  p  less  than  unity,  whereas  the 
stationary-phase  result  closely  approximates  the  exact  result  for  p  values  in 
excess  of  unity  as  stated  by  the  Equation  (72)  approximation. 

To  illustrate  the  degree-oi-aecuracy  of  the  filter-function  approximation 
as  it  affects  heterodyne  efficiency  calculations,  Figure  15  provides  a 
comparison  of  the  heterodyne  efficiency  calculated  by  the  exact  two-fold 
integral  analysis  [Equation  (48)]  with  the  corresponding  efficiency  calculated 
by  the  four-fold  integral  [Equations  (56)  and  (67)] ,  assuming  a  comparable 
Gaussian  source  LRCS  distribution.  The  Figure  15  data  demonstrate  that 
efficiency  results  can  be  computed  using  the  general  four-fold  integral  method 
which  closely  approximates  results  obtained  with  a  more  exact  (two-fold 
integral)  technique  of  limited  applicability. 

6.3*2  Atmospheric  Turbulence  Effects 

Figure  16  illustrates  the  predicted  effects  of  atmospheric  turbulence  on 
heterodyne  efficiency  for  (a)  a  glint  return  [Equations  (56)  and  (59)]  and  (b) 
the  return  from  a  small  (1  cm  radius)  spherical  target  [Equations  (56)  and 
(67)]  as  a  function  of  the  range  to  the  target.  The  Figure  16  data 
demonstrate  that  although  atmospheric  turbulence  decreases  the  heterodyne 
detection  efficiency,  especially  at  long  ranges  (~  10  km),  the  degradation  in 
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signal-to-noise  ratio  is  relatively  minor  compared  with  diffuse  target  speckle 
effects.  The  computed  results  plotted  in  Figure  16b  also  demonstrate  that  the 
heterodyne  efficiency  computed  for  the  small  sphere  (using  the  approximate 
four-fold  integral  method)  closely  approximates  the  glint  return  results, 
substantiating  the  validity  of  the  approximate  filter  function  for  turbulence 
effects  calculations. 


(a)  <b) 


OPOWMHJ 


Figure  13.  Laser  radar  filter  function  calculated  by  <a)  numerical  integration,  (b)  series 

expansion,  and  (c)  stationary  phase  approximation  as  a  function  of  the  unitiees 
variable  p  *  Ring’  and  the  parameter  (T «  Rm'^. 
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Figure  14.  User  radar  filtar  function  calculated  by  (a)  numerical  integration,  (b)  series 

expansion,  and  (c)  stationary  phase  approximation  as  a  function  of  the  unitless 
variable  p  ■  Rmp'  and  the  parameter  /f ■  Rm'^. 


0.2  0.3 

Target  radius  (m) 


GP03  0S38  IS 


Figure  15.  Calculated  laser  heterodyne  maximum  detection  efficiency  as  a  function  of 
the  radius  of  a  Schell-model  source.  Fixed  parameters  of  the  calculation  are 
A  =  P  =  2.5  cm,  f^.  =  ff  =  25  cm,  and  R  =  500  m. 


6.3.3  Target  Reflectivity  Calculations 

To  provide  a  data  file  of  representative  target  LRCS  distributions  for 
subsequent  studies  of  target  effects  on  heterodyne-detection  efficiencies,  a 
set  of  representative  (yet  hypothetical)  rough-surface  reflectivities  was 
defined  using  the  R0SSC0  laser  scattering  algorithm^,  and  various  projected 
LRCS  distributions  were  calculated  using  the  projection  equations  developed  in 
Section  5.4. 2»  Figure  17  depicts  the  backscatter  bidirectional  reflectance 
distribution  functions  that  are  predicted  by  the  ROSSCO  algorithm  for  metallic 
surfaces  having  a  small-scale  surface  roughness  described  by  an  exponential 
surface-height  power  spectral  density  function  (with  varying  correlation 
lengths)  and  a  Gaussian  surface-height  probability  density  function.  These 
functions  assumed  characteristics  of  the  small-scale  roughness  correspond  to 
measured  properties  of  many  machined  and  industrially  processed  metallic 
materials.  *  The  hypothetical  rough  surfaces  are  numbered  1  through  4  in 
increasing  order  of  smoothness.  Figure  18  provides  three-dimensional  graphs 
of  the  spherical-surface  LRCS  distribution  function,  Equation  (105),  for  the 
same  metallic  surfaces.  The  distributions  plotted  in  Figure  18  show  the 
transition  from  a  nearly  uniform  circular  LRCS  distribution  for  a  rough 
sphere,  Figure  18a,  to  the  point-like  distribution  for  a  very  smooth  sphere, 
Figure  18d. 
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Heterodyne  efficiency. 


(a) 

Glint  return 
2-D  integral  result 


Range  (km) 


(b) 

Glint  return  vs 
small  target  result 


Range  (km) 
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Figure  16.  Laser  heterodyne  maximum  detection  efficiency  computed  as  a  function  of 
range  for  (a)  a  single  target  glint  using  the  two-fold  integral  Schell-model 
method  and  (b)  a  1  cm  radius  target  (having  a  Gaussian  intensity  distribution) 
using  the  approximate  four-fold  integral  method.  Fixed  parameters  of  the 
calculation  are  A  =  P  =  2.5  cm,  f^  =  f,  =  25  cm,  and  LQ  =  5  cm. 
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Figure  17.  Laser  backscatter  bidirectional  reflectance  distribution  functions  for  hypo¬ 
thetical  metallic  rough  surface  materials  characterized  by  the  ROSSCO 
parameters  k a  =  3.0,  s  =  0.001,  m  =  1.0,  n  =  100.0,  a  »  0  and  correlation 
distances:  (a)  kK  *  1.0  (rough  surface  1),  k8  =  2.0  (rough  surface  2),  (b)  kC  =  5.0 
(rough  surface  3),  and  (c)  k?  =  10.0  (rough  surface  4). 


(a) 


(b) 


0.4 
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Figure  18.  Projected  laser  crou  section  per  unit  area  distributions  for  spherical  targets 

(0.2S  m  radius)  having  surface  characteristics  of  (a)  rough  surface  1,  (b)  rough 
surface  2,  (c)  rough  surface  3,  and  (d)  rough  surfaca  4. 

To  examine  the  effects  of  large-scale  surface  roughness,  two  different 
sets  of  roughness  parameters  were  postulated  that  result  in  the  Figure  19 
ROSSCO  predictions  of  the  backscatter  laser  bidirectional  distribution 
functions.  Large-scale  surface  undulations  produce  an  ensemble  of  surface 
glint  points  for  laser  reflectance  in  addition  to  the  dlffraction-like 
scattering  from  the  small-scale  roughness.  To  examine  the  limiting  effects  of 
large-scale  roughness,  a  surface  vas  defined  (rough  surface  A)  consisting 
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(a) 


(b) 
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Figure  19.  Laser  backscatter  bidirectional  distribution  functions  for  hypothetical  metallic 
rough  surface  materials  having  two  scales  of  roughness  characterized  by  the 
ROSSCO  parameters  (a)  ko  *  0.0,  kH  =  5.0,  k6  =  3.0,  s  =  0.3,  m  =  1.0,  n  =  100.0, 
a  =  0  {rough  surface  A),  and  (b)  ko  =  0.2,  kH  =  5.0,  k?  =  3.0  s  =  0.6,  m  =  1.0, 
n  =  100.0,  and  a  =  0  {rough  surface  B). 


solely  of  glint  points,  while  a  second  surface  was  defined  (rough  surface  B) 
having  a  much  larger  range  of  surface  undulations  combined  with  a  small-scale 
surface  roughness.  The  Figure  19a  data  show  the  rapid  decrease  in  the 
reflectivity  of  rough  surface  A  with  increasing  backscatter  angle  caused  by 
the  rapid  decrease  in  the  number  of  glint  points.  The  reflectivity  of  rough 
surface  B,  Figure  19b,  contrasts  with  that  of  rough  surface  A  because  the 
large  surface  undulations  and  small-scale  roughness  combine  to  produce  an 
almost  constant  reflectivity  distribution  (a  Lambertian  surface  has  a  constant 
value  of  p'  at  all  angles).  Figure  20  shows  the  LRCS  distribution  of 
spherical  targets  having  surface  features  parameterized  by  the  ROSSCO 
parameters  of  rough  surfaces  A  and  B.  The  principal  reflectance  from  the 
sphere  of  rough  surface  A,  Figure  20a,  arises  from  the  central  distribution  of 
glint  points,  whereas  significant  reflectance  from  the  rough  surface  B  sphere. 
Figure  20b,  occurs  at  locations  at  the  sphere  periphery. 
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(a) 


(b) 


Figurt  20.  Projected  lasar  cron  taction  par  unit  araa  distributions  for  spherical  targets 

(0.25  m  radius)  having  surface  characteristics  of  (a)  rough  surface  A  and  (b)  rough 
surfaces. 


To  illustrate  the  effects  of  surface  roughness  on  the  LRCS  distribution, 
of  conical  targets.  Figures  21  and  22  provide  calculated  LRCS/area  results, 
Equations  (107)— (1 14) ,  for  cones  of  rough  surfaces  A  and  B  at  various 
observation  angles.  A  comparison  between  the  data  plotted  in  Figures  21  and 
22  shows  the  expected  result  that  the  principal  reflectance  contribution  from 
cone  A  arises  from  an  area  near  the  cone  axis  (for  nonzero  observation  angles) 
as  a  result  of  the  abundance  of  glint  points  near  the  cone  axis,  whereas  the 
LRCS  of  cone  B  is  more  evenly  distributed  over  the  projected  outline  of  the 
cone  as  a  result  of  the  significant  laser  scattering  produced  by  rough  surface 
B  at  angles  near  grazing  incidence.  The  subsequent  section  demonstrates  the 
effect  of  varying  LRCS  distributions,  arising  from  surface  roughness 
differences,  on  calculated  heterodyne  detection  efficiencies. 
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Figure  21.  Projected  later  cross  section  per  unit  area  distributions  for  a  conical  target 

having  a  base  radius  of  0.25  m,  a  height  of  1.0  m,  and  a  surface  finish  character¬ 
ized  by  rough  surface  A.  Projection  direction  angles  are  (a!  0°,  (b)  20°,  <c)  40°, 
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Figure  22.  Projected  laser  cross  section  per  unit  area  distributions  for  a  conical  target 
having  a  base  radius  of  0.2S  m,  a  height  of  1 .0  m,  and  a  surface  finish 
characterized  by  rough  surface  B.  Projection  direction  angles  are  (a)  0°, 


6.3.4  Efficiency  Calculations 


A  variety  of  LRCS  distributions,  such  as  those  depicted  in  Figures  18, 

20,  21,  and  22,  were  used  to  determine  the  effects  on  the  calculated 
heterodyne  efficiency  of  (a)  target  roughness,  (b)  target  shape,  (c)  target 
orientation,  and  (d)  the  dimension  of  the  data  matrix  used  in  the  numerical 
integration  of  Equation  (67).  Calculated  values  of  the  heterodyne  efficiency 
for  spherical  targets  having  rough  surfaces  corresponding  to  rough  surfaces  1- 
4  are  tabulated  in  Table  1  for  various  dimensions  of  the  variable  grid  size 
employed  in  the  numerical  outer  integration  (dummy  variable  f+)  and  inner 
integration  (dummy  variable  of  Equation  (67).  The  Table  1  data  show  (1)  a 

monotonic  increase  in  the  calculated  heterodyne  efficiency  with  increasing 
smoothness  of  the  surface  and  (2)  relatively  little  effect  of  the  calculation 
grid  size  on  computed  results,  although  the  relative  error  resulting  from 
small  grid  sizes  increases  with  increasing  surface  roughness.  Both  of  the 
trends  noted  in  the  Table  1  data  are  expected  from  physical  and  mathematical 
considerations;  however,  it  is  noteworthy  that  (a)  solely  surf ace- roughness 
effects  can  account  for  predicted  heterodyne  efficiencies  differing  by  an 
order  of  magnitude  for  the  same  target  and  (b)  computational  grid  sizes  as 
small  as  5x5  can  yield  efficiency  results  that  are  in  relatively  good 
agreement  with  results  computed  from  grid  sizes  of  much  higher  dimensionality. 


TABLE  1.  ROUGH  SURFACE  AND  GRiD  SIZE  EFFECTS  ON  CALCULATED 
HETERODYNE  EFFICIENCY 


Fixed  parameters  of  the  calculation  are  A  =  P  =  2.5  cm,  f^, 
R  =  500  m,  Cn2  =  0,  sphere  radius  =  0.25  in 


kV 

Rough 

surface 

Grid  size 

5x5 

(5x51 

9x9 

(9x9) 

11  x  11 
(11  x  11) 

5x11 
(11  x  11) 

5x  11 
(5  x  11) 

1 

1 

0.13 

0.10 

0.11 

0.10 

0.11 

2 

2 

0.37 

0.27 

0.27 

0.28 

0.30 

5 

3 

0.95 

0.82 

0  74 

0.83 

0.88 

10 

4 

1  00 

1  00 

0.95 

0.97 

1.00 

=  fa  25  cm, 


— —  Innt '  integration 
-* —  Outer  integration 


Note  Rounded  to  two  significant  figures 
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Calculated  spherical-target  heterodyne  efficiencies  are  tabulated  in 
Table  2  for  various  ranges  assuming  rough-surface  textures  corresponding  to 
rough  surfaces  A  and  B  and  a  centrally  located  spherical  glint  point.  The 
Table  2  data  demonstrate  the  monotonic  increase  in  the  heterodyne  efficiency 
with  increasing  range  to  the  target  and  again  illuscrate  the  reduced 
heterodyne  efficiency  resulting  from  rough  surfaces.  It  is  apparent  from  the 
Table  2  data  that  glint  points  are  relatively  more  effective  in  enhancing  the 
heterodyne  efficiency  of  rough  targets  vis-a-vis  smooth  targets. 

Calculated  conical  target  heterodyne  efficiencies  are  plotted  in  Figure 
23  as  a  function  of  range  for  rough  surfaces  1-4  and  in  Figures  24-25  as  a 
function  of  aspect  angle  for  rough  surfaces  A  and  B.  The  cone  efficiency  data 
plotted  in  Figures  23-25  substantiate  the  conclusions  drawn  from  the  spherical 
target  calculations  regarding  the  deleterious  effects  on  the  heterodyne 
efficiency  of  roughened  surfaces  and  decreased  range.  Additionally,  however, 
the  data  of  Figures  23-25  demonstrate  that  the  aspect  angle  of  three- 
dimensional  bodies  critically  influences  the  heterodyne  efficiency  because  of 
the  variation  in  the  projected  area  with  asoect  angle. 

As  a  final  demonstration  of  the  effects  of  the  grid  size  used  in  the 
numerical  calculation  of  the  heterodyne  efficiency.  Table  3  provides  a 
comparison  of  the  range-dependent  heterodyne  efficiencies  calculated  for  a  60° 
aspect  angle  cone  of  rough  surface  B  for  two  different  grid  sizes  (11x11  and 
21x21).  The  relatively  minor  discrepancies  observed  in  the  Table  3  data  for 
the  different  grid  size-calculations  indicate  that  sufficiently  precise 
efficiency  computations  can  be  performed  economically  for  relatively  diffuse 
scattering  targets.  Because  the  major  part  of  the  computation  time  of  the 
heterodyne  efficiency  is  devoted  to  the  four-fold  integral,  Equation  (67),  it 
should  be  expected  that  the  increased  running  time  of  the  21x21  grid  size  of 
Table  2  over  the  lixll  grid  is  given  by  the  factor 

(21/1 1 )4  -  13.28, 

which  is  in  good  agreement  with  the  actual  factor  of  a  12.6  increase  in 
running  time  indicated  in  Table  3.  Thus  it  is  expected  that  100  s  of 
computation  on  the  CDC  6500  (the  machine  used  for  the  calculations  of  this 
study)  are  sufficient  to  perform  a  heterodyne  efficiency  computation  on  a  grid 
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size  of  38x38  (both  inner  and  outer  integrations)  vrtiich  is  sufficient  to 
estimate  the  heterodyne  detection  efficiency  resulting  from  a  rather  complex 
geometrical  shape* 


TABLE  2.  SPHERICAL  TARGET  HETERODYNE  EFFICIENCIES 

•  Receiver  parameters:  A  =  P  =  25  cm,  fg  =  fr  =  25  cm,  L0  =  5  cm 

•  Target/encounter  parameters:  radius  =  25  cm,  a^  =  1 ,  Cn^  =  10'^4  m'2/3 


Range 

(km) 

Heterodyne  efficiency 
incoherent/total 

Rough 
surface  A 

Rough 
surface  B 

0.5 

0.59/0.78 

0.18/0.53 

1.0 

0.87/0.93 

0.64/0.74 

10.0 

0.94/0.96 

n.93/0.95 

CiP03  0«  39-27 


TABLE  3.  GRID  SIZE  EFFECT  ON  CALCULATED  LASER  DETECTION 
HETERODYNE  EFFICIENCIES  FOR  A  CONICAL  TARGET 


F'xed  parameters  of  the  calculation  are  A  =  P  “  2.5  cm,  fg  =  fr  =  25  cm,  L0  =  5  cm, 
=  10"^  m'2  ,  and  60°  projection  angle  for  a  cone  having  a  surface  finish 
characterized  by  rough  surface  B. 


Grid  size 

LRCS 

Heterodyne  efficiency 

Run 

time 

(s) 

Range  (km) 

0.5 

1.0 

5.0 

10.0 

11  x  11 

0.291 

0.0364 

0.0719 

0.695 

0.915 

0.71 

21  x  21 

0.285 

0.0317 

0.0650 

0.688 

0.907 

8.96 

%  Difference 

2 

15 

11 

1 

1 

- 
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Heterodyne  efficiency. 


Figure  23,  Laser  heterodyne  maximum  detection  efficiency  values  computed  as  a  function 
of  range  for  coi  ical  targets  having  a  base  radius  of  0.25  m,  a  height  of  1.0  m, 
and  the  annotated  surface  finish  characteristics  and  projection  angles.  Fixed 
parameters  of  the  calculation  are  A  =  P  =  2.5  cm,  ?<,  =  f  =  25  cm,  L  =  5  cm, 
and  Cn2  =  10-K  m'2/3. 
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para  neters  of  the  calculation  are  A  =  P  =  2-5  cm,  fv  ,  •  a 

r  2  -  io-14nV2/3  and  range  values  are  annotated. 
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Figure  25.  Laser  heterodyne  detection  efficiency  values  computed  as  a  function  of  view¬ 
ing  (projection)  angle  for  a  conical  target  having  a  base  radius  of  0.25  m,  a 
height  of  1.0  m,  and  a  surface  finish  characterized  by  rough  surface  B.  Fixed 
parameters  of  the  calculation  are  A  =  P  -  2.5  cm,  =  fr  =  25  cm,  L0  =  5  cm, 
=  10'^  m '2/3,  and  range  (R)  values  are  annotated. 
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APPENDIX  A:  STATIONARY-PHASE  EVALUATION  OF  THE  FILTER  FUNCTION 


The  convolution  theorem  of  Fourier  transformation  can  be  used  to  write 
the  filter  function  [Equation  (63)]  in  the  form 
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upon  introducing  the  variable  change 
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the  integral  of  Equation  (A2)  can  be  approximated  by  the  method  of  double¬ 
stationary  phase^  because  the  complex  exponential  term  in  the  integrand 

dominates.  The  stationary  phase  point,  p  ,  is  located  at 
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The  double-stationary-phase  integral  formula  yields 


providing  Equation  (A4)  holds. 
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APPENDIX  Bs  SERIES  EXPANSION  EVALUATION  OF  THE  FILTER  FUNCTION 


Because  the  Heaviside  function  can  be  written  in  integral  form  as 
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the  filter  function  [Equation  (63)]  can  be  written  as 
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after  integrating  over  the  variable  t  •  A  contour  in  the  complex  t  plane 
bounded  by  the  real  axis,  and  an  infinite  arc  in  the  upper  half  plane  can  ba 
used  to  perform  a  contour  integration  of  Equation  (B3).  This  contour  encloses 
a  simple  pole  at  t  *  ie,  and  a  combination  pole  and  essential  singularity 
at  T  *  iB/2.  The  contribution  of  the  pole  and  essential  singularity 
combination  is  evaluated  by  noting  that 
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for  any  analytic  function  f(z).  Equation  (B4)  follows  from  the  Cauchy 
integral  theorem  upon  expanding  the  exponential  and  integrating  term  by  term 
using  the  contour  integral  method.  Also  noting  that 
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It  is  clear  that 
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providing 
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Thus,  utilizing  the  Equation  (B6)  aproximation  to  evaluate  the  combination 
pole  and  essential  singularity  contribution  to  the  contour  integration  of 
Equation  (B3),  the  approximate  result  is 
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providing 
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